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ABSTRACT 


We  study  the  equation 


pu  =  3n 
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where  u(x,t)  is  a  real  valued  function  of  x  e  [-1,1]  and  t  f  R,  with 
the  boundary  condition 
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at  x  =  +1.  This  equation  is  derived  as  a  model  for  the  elongation  of  thin 
filaments  of  polymeric  liquids,  u  denoting  the  position  of  a  fluid  particle 
in  space,  a  the  memory  kernel,  and  f  the  force  acting  on  the  ends  of  the 
filament.  We  study  the  evolution  of  u,  assuming  the  initial  condition 
u(x,t  =  -°°)  =  x.  It  is  shown  that  under  appropriate  conditions  on  a  and 
f  the  boundary  condition  can  be  uniquely  resolved  with  respect  to  u^.  The 
full  problem  is  transformed  in  such  a  way  that  it  is  approachable  by  the 
Sobolevskii  theory  of  quasilinear  parabolic  equations.  This  yields  the 
existence  of  solutions  to  the  initial  value  problem  on  sufficiently  small 
time  intervals.  Moreover,  we  show  that  if  f (t)  converges  to  zero  exponen¬ 
tially  as  t  -+■  +°°  and  is  small  in  an  appropriate  norm,  there  exists  a  solu¬ 
tion  globally  in  time,  which  approaches  a  stationary  limit  as  t  -*•  +°°. 
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SIGNIFICANCE  AND  EXPLANATION 


f  *  « -  ' '  -  j  f 

The  equation  under  study  -  stated  in  the  abstract  and  derived  from  physical 
principles  in  this  paper  -  describes  the  elongation  of  a  filament  of  a 
polymeric  liquid  subjected  to  a  force  f  at  both  ends.  Hie  liquid  is 
assumed  to  satisfy  certain  accepted  k rubberlike  liquid^  constitutive  rela¬ 
tions,  and  the  filament  is  assumed  to  be  thin,  which  permits  a  reduction  of 
the  problem  to  one  space  dimension.  The  unknown  variable  u  denotes  the 
position  of  a  fluid  particle  at  time  t,  which  was  at  position  x  at 
t  =  -°°,  i.e.,  before  the  deformation  started,  we  have  u(x,-“>)  =  x.  In 

this  paper  the  equation  under  study  is  transformed  in  such  a  way  that  it 
fits  into  the  framework  of  the  general  mathematical  theory  for  "quasilinear 
parabolic  equations'".  This  makes  it  possible  to  prove  that  for  any  given 
"initial  condition"  a  solution  exists  at  least  on  a  certain  time  interval. 
(It  is  a  part  of  the  analysis  to  discover  what  is  an  appropriate  meaning 
of  "initial  condition"  to  be  associated  with  the  problem  under  study). 
Moreover,  we  shall  prove  that  for  forces  f  (t) ,  which  approach  zero 
exponentially  for  t  *  +®  and  are  small  in  a  suitable  sense,  there  is  a 
solution  for  all  times  t,  -«®  <  t  <  +*  ,  and  this  solution  approaches  a 
stationary  limit  as  t  -*■  +°°. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  QUASILINEAR  PARABOLIC  EQUATION  DESCRIBING  THE 
ELONGATION  OF  THIN  FILAMENTS  OF  POLYMERIC  LIQUIDS 


M.  Renardy 

0.  INTRODUCTION. 

We  study  the  following  problem  occurring  in  polymer  processing:  A  thin  filament 
of  a  viscoelastic  liquid  is  subjected  to  a  force  f  acting  on  its  ends  as  shown  in 
the  diagram: 

f  4— — . .  . . .  ~  =3->  f 

We  investigate  the  temporal  evolution  of  the  displacement.  The  equations  that  our 
analysis  is  based  on  involve  the  “rubberlike  liquid"  constitutive  assumption  for  the 
stress-strain  law  [31  and  certain  approximations  based  on  the  thinness  of  the  filament, 
which  allow  the  reduction  to  a  spatially  one-dimensional  problem.  Using  these  assump¬ 
tions,  we  shall  derive  the  following  equation 


pu  =  3n 


3x3t 


.  ,  t  fu  (t)  u  (s)’1 

/  *'1  -  -tr  - -fd 

X  '■u  (s)  u  (t)' 


ds 
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X  X 

where  u(x,t)  is  a  real  valued  function  of  x  e  [-1,1]  and  t  t  3R.  As  usual,  a  sub¬ 
script  x  denotes  partial  differentiation  w.r.  to  x  and  "dot"  denotes  partial 
differentiation  w.r.  to  t.  The  arguments  (x,t)  are  suppressed  unless  needed  for 
proper  understanding.  (0.1)  is  supplemented  by  the  nonlinear  Neumann  boundary  condition 


(“)  +  /  a(t 


s) 


u  (t)  u  (s) 
x  x 


uds)  U2(t) 
X  X 


ds  =  f(t) 


(0.2) 


at  x  =  +1. 


In  these  equations  u(x,t)  denotes  the  position  at  time  t  of  a  fluid  particle, 
which  is  at  the  position  x  in  a  certain  reference  state.  This  reference  state  will 
be  identified  with  the  state  of  the  fluid  at  t  =  -<»,  i.e.  we  have  u(x,t  =  -°°)  =  x. 

0  denotes  the  density  of  the  fluid,  n  the  viscosity,  and  f  the  force  acting  on  the 
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ends  of  the  filament.  The  memory  kernel  a  :  [0,°°)  -*■  IR  will  be  assumed  to  have  the 
following  properties,  which  we  shall  refer  to  as  assumptions  (a): 

(i)  a  has  the  representation 

a  (t)  -  /  e  Xtdp(X)  (0.3) 

where  p  is  a  complex  valued  Borel  measure  on  the  complex  plane  <C  such  that 
1  f  LX(u)  (i.e.  1  is  integrable  w.r.  to  the  total  variation  of  p) ,  and  supp  p  is 

contained  in  {X  £  <E  |  — _<  arg  X  ,  |x|  >_  t}  for  some  ^  <  -j  and  c  >  0.  Since  a 

is  real,  we  may  and  will  assume  that  dp(X)  -  dp(X). 

(ii)  aft)  ^  0  for  t  £  (0,»). 

(iii)  a  is  monotonely  decreasing. 

Note  that  (i)  implies  in  particular  that  a  is  continuous  and  that  |a|  can  be 
estimated  by  a  decaying  exponential.  The  motivation  for  assumptions  (a)  will  become 
apparent  later  in  the  paper.  The  sectorial  condition  for  supp  p  is  needed  to  make  the 
problem  fit  into  the  theory  of  parabolic  equations,  (ii)  and  (iii)  will  have  important 
implications  for  the  spectra  of  certain  linear  operators.  In  physical  theories  derived 
from  "molecular  network”  or  "bead-spring"  models  (see  14)  and  the  references  in  13] , 
ch.  6)  a  turns  out  to  be  a  finite  sum  of  decaying  exponentials.  This  is  clearly  a 

special  case  of  assumptions  (a),  p  in  this  case  being  a  finite  sum  of  Dirac  measures 

located  on  the  real  axis. 

The  boundary  condition  (0.2)  agrees  precisely  with  the  equation  describing  the 
evolution  of  the  length  of  the  filament  when  inertial  forces  are  neglected.  This 
problem  has  been  discussed  previously  by  lodge,  McLeod  and  Nohel  in  [5j  and  by  the 
author  in  [7] .  Lodge,  McLeod  and  Nohel  consider  the  solution  as  known  for  t  <  0  and 
assume  it  is  nondecreasing.  They  then  assume  f  »  0  for  t  >  0  and  study  existence, 
asymptotic  behaviour  and  various  monotonicity  properties  of  solutions.  In  [7]  the 
force  f  is  a  giver,  continuous  function  IR  IR.  It  is  assumed  that  f  converges  to 
zero  exponentially  as  t  *  and  that  either  f  converges  to  0  exponentially  as 

t  -*  +*»  and  is  small  in  a  suitable  norm,  or  the  size  of  f  is  arbitrary,  but  f(t) 


vanishes  identically  for  t  greater  than  some  finite  t^.  (In  the  latter  case  we  need 

the  additional  assumption  that  supp  u  is  on  the  real  axis)  in  fact  in  [7]  we  assumed 

that  u  was  a  finite  sum  of  i'irac  measures  on  the  real  axis,  but  the  same  ideas  can 

be  applied  to  the  more  general  situation  as  we  demonstrate  below),  in  both  cases  it  is 

proved  that,  given  the  initial  condition  (t  ■  -“)  »  1,  a  unique  positive  solution 

exists  for  all  times  t,  and  moreover  lim  u  (t)  exists  and  is  strictly  positive. 

t-++«>  * 

Whereas  the  arguments  in  (5]  rely  mainly  on  monotonicity  properties,  the  main  tools  in 
(7]  are  the  implicit  function  theorem  and  the  use  of  Liapunov  functions. 

The  present  paper  will  be  arranged  as  follows:  In  51  we  explain  the  basic  physical 
laws  and  the  approximations  leading  to  (0.1),  (0.2).  We  start  from  the  basic  laws  of 
continuum  mechanics,  using  the  "rubberlike  liquid"  constitutive  relation.  The  equation 
of  motion  in  the  interior  of  the  fi lament  and  the  boundary  conditions  on  the  lateral 
surface  are  then  solved  formally  by  a  power  series  expansion  with  respect  to  a  "thinness 
parameter"  in  an  analogous  manner  as  was  done  in  the  theory  of  thin  elastic  rods  (6) . 

The  first  order  terms  in  this  expansion  lead  to  (0.1).  The  formal  expansion  does  not 
in  general  fit  given  boundary  conditions  at  the  ends  of  the  filament,  and  one  is  con¬ 
fronted  with  a  "boundary  layer"  problem.  Since  we  are  only  interested  in  a  first  order 
approximation,  we  shall  not  deal  with  this  situation  here.  Instead,  we  consider  the 
balance  of  force,  taking  into  account  only  terms  not  involving  the  small  parameter. 

This  leads  to  (0.2).  §2  summarizes  the  results  of  [7]  concerning  (0.2)  as  explained 

above,  taking  into  account  the  modifications  required  by  the  more  general  assumptions 
on  a.  As  a  result,  we  may  subsequently  consider  ux  as  being  given  on  the  boundary. 

In  53  and  54  we  finally  deal  with  the  full  problem  (0.1),  (0.2).  Using  (0.3),  this 
problem  is  transformed  in  such  a  way  that  it  fits  into  the  abstract  theory  of  quasi- 
linear  parabolic  equations  introduced  by  Sobolevskii  S2),  IS).  An  "initial  condition" 
in  th-~  evolution  problem  so  defined  will  not  necessarily  involve  the  whole  history  of 
u,  but  only  certain  of  its  moments,  the  choice  of  which  depends  on  the  support  of  y. 

In  53  we  shall  explain  this  transformation  and  as  a  consequence  of  the  Sobolevskii 


theory  obtain  the  existence  and  uniqueness  of  solutions  to  the  initial  value  problem 
locally  in  time.  §4  deals  with  the  case  where  f  converges  to  zero  exponentially  as 
t  ■+  +~  and  is  small.  We  assume  the  filament  is  undeformed  (u  »  x)  at  t  =  It 


will  be  shown  that  a  solution  of  the  full  problem  exists  globally  in  time,  which 
approaches  a  stationary  limit  as  t  ■*  +». 


1.  DERIVATION  OF  THE  BASIC  EQUATIONS. 

We  assign  to  each  point  in  the  fluid  two  different  sets  of  coordinates:  By 
12  3 

(S  ,  C  ,  G  )  we  denote  "body  coordinates",  i.e.  coordinates  labelling  a  specific  particle 

in  the  fluid.  These  coordinates  can  be  identified  with  the  position  of  the  particle  in 

space,  when  the  fluid  is  in  a  certain  "reference  state".  (It  will  later  be  convenient 

to  take  as  a  reference  state  the  state  of  the  fluid  at  time  t  =  -<») .  On  the  other 

hand  (y  ,y  ,y  )  will  denote  coordinates  labelling  a  point  in  space.  We  are  interested 

in  finding  trajectories  of  fluid  particles,  i.e.  a  functional  dependence  y  (G  ,G  ,G  ,t). 

In  our  exposition  of  the  equations  describing  this  functional  dependence  we  follow 

Dodge  (3].  (For  a  summary  of  the  relevant  equations,  see  p.  206-207). 

12  3 

To  each  point  (S  ,?  ,z,  )  there  is  assigned  a  body  metric  tensor  y  and  a  body 
stress  tensor  it.  y  is  defined  by  the  relation 

V  =  i4  .  hi 

3C1 

ir  is  related  to  y  by  a  constitutive  law,  which  expresses  the  specific  properties  of 
the  material.  We  use  the  "rubberlike  liquid"  constitutive  relation  ([3],  p.  143): 

tt1^  +  py1^  =  -n  ■  +  /  a  (t  -  s) y1^  (s)ds 


where  y*^  denote  the  components  of  y  x,  i.e.  =  <5^;  n  is  a  positive 

material  constant  called  the  viscosity,  and  a  is  a  given  function,  which  will  always 
be  assumed  to  satisfy  assumptions  (a)  stated  in  the  introduction,  p  is  an  unknown 
variable  having  the  physical  significance  of  a  pressure.  The  introduction  of  this 
variable  is  necessary,  since  we  assume  the  fluid  is  incompressible 

det  y  -  1  (1.1) 

k  • 

The  evolution  of  y  is  determined  by  Newton's  law,  which  for  a  Cartesian  space 
coordinate  system  takes  the  form 

k 

L. 

3r.X  l3G" 


0y*  =  *4  U-  *si  +  r1  ,rs| 

3,1  V  rs  > 


(1.2) 


-5- 


The  Trs  denote  the  Riemann-Christof f el  symbols  associated  with  the  metric  tensor  y. 


r1 

rs  2  \  a 


f3Y-‘r  3Yjs  3y- 
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Equations  (1.1)  and  (1.2)  have  to  be  supplemented  by  boundary  conditions  referring  to 


f) 


either  the  displacement  or  the  stresses  on  the  boundary  of  the  liquid,  we  shall  here 

]( 

deal  with  stress  conditions.  Let  Y  denote  the  components  of  surface  traction 

l 

referred  to  space  coordinates.  Then  the  boundary  conditions  on  a  surface  5  «  const, 

are  given  by 


*il  *4  (v“)'1/J  =  Yk  (1.3) 

351 

In  the  problem  of  the  elongated  filament  the  surface  traction  on  the  lateral  surface 
is  zero,  whereas  at  the  ends  there  is  a  longitudinal  surface  traction  equal  to  f 

divided  by  the  cross-sectional  area  of  the  filament.  For  convenience,  we  let  5*  and 

1  2  3 

y  resp.  denote  the  coordinate  is  the  direction  of  the  filament  and  5  ,5  and 

2  3 

y  , y  resp.  the  transversal  coordinates.  It  is  assumed  that  in  the  undeformed  refer¬ 
ence  state  (i.e.  at  t  =  -»)  the  filament  is  cylindrical  and  axi -symmetric,  i.e.  in 

this  state  we  have  y^  =  ,  where  5*  (by  appropriate  normalization  of  length  scale) 

/  2  2  3  2 

ranges  from  -1  to  1,  and  r  =  /(5  )  +  U  )  ranges  from  0  to  5,  the  radius 

of  the  filament.  Then,  for  small  6,  equations  (1.1)  and  (1.2)  and  the  boundary 

conditions  on  the  lateral  surface  can  formally  be  solved  by  a  series  expansion  in 

powers  of  r  and  S.  This  expansion  is  analogous  to  that  used  by  Nariboli  [6]  for 

the  problem  of  longitudinal  elastic  waves  in  a  thin  rod. 

2  3 

We  put  t  =  7—,  r  =  t,  so  that  the  lateral  surface  now  corresponds  to 

0  0  0 

r  =  1.  We  make  the  following  ansatz 


-6- 


2,1  xp2  . r  .2\>  1  ~2 

y  (C  ,5c  ,5c  >  -  6?  J.  5  Q  (C  ,r  ) 
vO 


y3(c\6C2,«C3>  -  «C3  [  62vo  U1^2) 

v-0 


P(c\«c2.«c3>  =  £  «2vr  u1,^2) 

v»0  v 

A  2 

where  P^Q^.R^  are  polynomials  of  vth  degree  in  r  .  This  ansatz  is  inserted 
into  (1.1),  (1.2)  and  the  lateral  boundary  conditions,  which  are  supposed  to  be 
satisfied  for  all  values  of  6.  Formally  this  yields  an  infinite  set  of  equations 
for  the  coefficients  of  PV«QV  and  Ry.  We  are  only  interested  in  deriving  an  equa¬ 
tion  for  the  first  term  P^fc^),  and  we  shall  in  the  following  only  carry  out  the 
series  expansion  as  far  as  needed  for  this  purpose. 

when  terms  up  to  0(6)  are  taken  into  account,  we  find  for  the  metric  tensor 


(  r^oi2 

3C1 


si2-* 


6?  ■* 


-2  *3 

6t  •*  6c  •* 


where 


Z!o  !i:-l 

.1 


3  Q 
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3;  3r  3c 

Putting  6=0,  we  find  from  (1.1) 


3V2 


3C 


1)  C0  =  1 


(1.4) 


Next  we  consider  the  boundary  conditions  on  the  lateral  surface.  At  a  boundary  point 

-3  21  22  23 

where  r,  =  0  these  yield  the  equations  n  =  tt  =  n  =0.  (Because  of  the  radial 

symmetry  it  suffices  to  consider  these  boundary  points;  if  the  traction  vanishes  there. 


-7- 


23 

it  does  so  everywhere  else.)  it  vanishes  identically  as  a  result  of  the  radial 
22 

symmetry.  For  it  we  obtain  the  following  terms  of  the  order  0(1) 

n22  =  -RgQ~2  -  n  (Q02)  +  /  a(t  -  s)Q02(s)ds  =  0  (1.5) 

—  00 

Finally  we  have 


21 


-PY 


21 


3  .  21 

at  lY  } 


+  /  alt  -  s)y21(s)ds  = 


(1.6) 


0 

All  solutions  we  are  going  to  consider  shall  satisfy  lim  pit)  =  /  a(-s)ds  and 

+•-►—00  _00 

21 

lim  Y  (t)  =  0,  the  convergence  being  exponential.  If  (a)  holds,  it  is  then  not 
t~>—  «© 

difficult  to  prove  that  the  only  solution  to  (1.6)  satisfying  the  specified  conditions 
21 

is  y  =0.  In  the  first  order  in  5  this  yields  =  0.  The  law  of  motion  (1.2) 
now  yields  the  following  equation  for  Pq|y^  = 


In  order  to  simplify  notation,  we  shall  henceforth  write  u  for  PQ  and  x  for  C*. 
The  last  equation  now  yields  (0.1)  after  a  few  manipulations,  when  (1.5)  and  (1.4)  are 
used  to  express  in  terms  of  u. 

Finally  we  have  to  specify  boundary  conditions  at  the  ends  of  the  filament.  As 
noted  in  [6],  the  asymptotic  expansion  which  we  used  for  the  interior  problem  generally 
fails  near  the  ends,  and  a  "boundary  layer”  has  to  be  taken  into  account.  The  boundary 
layer  is  discussed  in  a  forthcoming  paper  by  Reiss,  which  is  referenced  in  [11,  but  not 
available  yet.  We  are  here  only  concerned  with  a  first  order  approximation,  and  we 


shall  ignore  boundary  layer  effects.  Instead,  we  take  care  of  the  force  balance  in 
the  zeroth  order  with  respect  to  6.  Namely,  if  one  formally  inserts  our  expansion 


into  the  boundary  conditions  at  the  ends,  it  is  seen  that  ail  traction  components 
transversal  to  the  direction  of  the  filament  are  0(6).  The  longitudinal  traction 
component  gives  the  following  terms  of  order  0(1) : 


11  3P0  ,  11 . -1/2 
t  =  n  — r  (y  ) 

H 


ii  f3po)2 


U  t  fU  (t)  U  (S)i 

3n  -*■  +  /  a (t  -  s)  -f - * 

u  *  2  u 

X  -®  VU  (s)  X 


it)  J 


ds 


Since  the  cross-sectional  area  of  the  filament  is  in  first  approximation  equal  to  u 


-1 


we  shall  require  that  r  =  f*u  .  This  yields  (0.2). 


2.  THE  BOUNDARY  PROBLEM. 


i 


I 


I 


In  this  section  we  consider  the  problem  of  solving  (0.2)  for  u^,  when  f  is 
given.  The  results  we  present  slightly  generalize  those  of  (7) ,  allowing  for  the  more 
general  class  of  kernels  a  satisfying . assumptions  (a).  In  order  to  simplify  notation, 
we  write  y  for  u^.  Instead  of  (0.2)  we  study  the  slightly  more  general  problem 


3ny  +  /  a(t  -  s)  Of  ***  -  y(s)]ds  »  f(t)ya 

ly  (s)  > 


(2.1) 


where  0  <  a  <  3.  (0.2)  corresponds  to  a  »  2,  and,  as  explained  in  (7),  the  case 

a  =  is  also  physically  interesting,  namely,  it  describes  the  deformation  of  a  sheet 
of  the  polymer,  when  inertia  are  neglected: 


/"  i 


y 


-f 


/ 


-f 


We  put  g  (X)  =  /  e  *  ^  S'y  2(s)ds,  h(X)  =  /  e  S'y(s)ds, 

••00  —00 

Y  (X)  =  g(X)y2,  6 (X)  =  h(X)y_1 

Then  (2.1)  is  equivalent  to  either  of  the  following  systems 

3ny  =  /  (h  (X)  -  g  ( X ) y 3 ) dy  ( X )  +  fya 
g  (X)  =  -Xg(X)  +  y  2  (2.2) 

h(\)  =  -Xh(X)  +  y 

3ny  =  y  /  ( A  ( X )  -  y(X))du(X)  +  fya 

Y  (  X )  =  -Xy  (X)  +■  1  -  ~  Y(X)  /  (Y  ( A)  -  6(X))du(X)  +  y  (\)  fy*'1  (2.3) 

3n  3n 

A(X)  =  -  XA (X)  +  1  +  ~  MX)  /  (y  ( A )  -  5  { X) )  du  ( X)  -  i  S(X)fya_1 
3n  3n 

-lo- 


Both  forms  will  be  used  in  the  following.  Equations  (2.2)  or  (2.3)  will  be  regarded 

s  2  s 

as  evolution  problems  in  the  space  X  »  m  x  (L  (u))  (1  <_  s  <  “) .  Here  L  (u)  denotes 

the  space  of  all  (equivalence  classes  of)  functions  g  :  £  -*  C  such  that  g  (z)  ■  g  (z) 

and  |g | S  is  integrable  w.r.  to  the  total  variation  of  u.  Clearly,  the  right  side  of 
(2.2)  or  (2.3)  is  the  sum  of  an  analytic  generator  and  a  smooth  nonlinear  term. 

A  trivial  solution  for  f  »  0  is  given  by  y  «  1,  g(X)  «  h ( A )  «  -y(X)  »  5 ( X)  »  i, 
and  we  are  interested  in  solutions  converging  to  this  trivial  solution  as  t  *  As 

a  first  step  we  investigate  the  spectral  properties  of  the  linearization  of  (2.2)  (of 
course  (2.3)  gives  the  same  result)  at  this  point,  i.e.  we  study  the  inhomogeneous 
linear  equation 

3n6y  -  /  (h(X)  -  g(X))dy(X)  +  3y  /  j  du(X)  -  3nv> 

Bg(X)  +  Xg  (X)  +  2y  -  i^tX)  (2.4) 

Bh(X)  +  Xh(X)  -  y  -  ^(X) 

If  -B  is  not  in  the  support  of  u,  the  last  two  equations  can  be  resolved  with 
respect  to  g(X)  and  h(X).  This  inserted  into  the  first  equation  of  (2.4)  yields 

3n8y  -  3y  /  •— -g-  dy(X)  +  3y  /  i  dy(X) 

*  (X)  -  U-  (X) 

=  3rw>  +  /  — - -  du  (X) 

Hence  the  resolvent  exists  at  8,  iff  -B  is  not  in  the  support  of  u  and 


>  (3)  :=  3nP  -  3  /  (^-g-  -  i)du(X)  /  0 


Clearly,  p (8)  vanishes  for  8=0.  Namely,  we  have 


0(61  ”  36(11  +  /  uT7TTdll(3,) 


Using  the  relationship  between  u  and  the  kernel  a,  we  find 


/Mr77)d"(A>  =  / 


r 


For  S  J*  0,  the  real  part  of  this  expression  is  given  by 

CO 

/  Jgj2  ^  6(1  -  e  ^cos  (t  Im  g) )  +  I»  g  e_t  Re  ®sin  (t  im  6)  )dt 

If  Be  8^0,  condition  (a)  (ii)  implies  that  the  first  contribution  is  positive,  and 
condition  (a)  (iii)  implies  that  the  second  contribution  is  positive,  too.  Hence 
Re  /  g)du(A)  ^  o,  whence  certainly  p(B)  ■+  0. 

a  2 

For  easier  reference,  let  us  put  Y  »  (y,g,h)  e  R  *  (L  (u) )  in  (2.2)  and  write 
(2.2)  in  the  form 

y  »  MV  -  Vg)  +  N (Y  -  YQ,f)  (2.5) 

where  L  denotes  the  linearization  of  the  right  side  at  the  trivial  solution 
VQ  »  (l,  y, -j) .  Analogously,  we  put  y*  «  (y . Y , <5)  and  write  (2.3)  in  the  form 

Y*  «  if  (Y*  -  Yg)  +  N1  (Y«  -  Yg,f)  (2.6) 

we  have  just  proved 
Proposition  2.1; 

The  spectrum  of  L  (or  L' )  consists  of  the  algebraically  simple  eigenvalue  0 
(geometric  simplicity  is  immediate,  and  algebraic  simplicity  follows  from  the  fact  that 
the  resolvent  has  a  first  order  pole)  and  a  remainder  contained  in  the  left  half  plane. 
Moreover,  the  restriction  of  L  to  the  range  of  L  generates  an  analytic  semigroup 
of  negative  type. 

Before  we  can  state  our  theorems,  we  must  first  define  some  spaces  of  functions. 
Definition  2.2: 

bet  Z  be  a  Banach  space  and  a  a  positive  real  number.  Then 

X°(Z)  ;»  (v  e  Cn(ffi,  Z)  j  lim  e°  I t  i  ||v (k>  (t)  I!  -  0  for  k  =  0,1, . . .  ,n  ; 

"  tH+» 

(k) 

v  denoting  the  kth  derivative} 


l 
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function  argument  as  above,  but  ratber  than  considering  functions  on  all  of  ®,  we 
consider  only  function  on  (-“.t^l.  From  this  we  see  that  Y  is  unique  in  the  class 
of  all  functions  that  approach  zero  exponentially  as  t  -*  Finally,  if  we  assume 

Y  converges  to  zero  at  all,  it  can  be  seen  from  the  last  two  equations  of  (2.3)  that 
v  -  y  and  5  -  i  converge  to  zero  exponentially,  because  if  only  these  two  equations 
are  considered,  the  zero  eigenvalue  in  the  linearization  does  not  occur.  From  the  first 
equation  of  (2.3)  we  find  that  y  converges  to  zero  exponentially,  and  hence  the  con¬ 
verge  of  y  to  its  limit  has  to  be  exponential,  too. 

If  further  restrictions  are  made  on  u,  *  global  result  can  be  proved  that  does 
not  rely  on  the  smallness  of  f. 

T.eorem  2.4; 

In  addition  to  (a),  assume  supp  u  is  contained  in  the  real  axis  and  u  is  posi¬ 
tive  real.  Let  f  :  IR  -»  TR  be  continuous  and  such  that  lim  e  otf(t)  •  0  for  some 

t—  •> 

■>  ■>  0  and  f(t)  *  0  for  t  >  For  any  such  f,  equation  (2.3)  has  a  unique 

solution  satisfying  lim  Y* (t)  «  yq.  This  solution  exists  globally  in  time,  moreover, 

lim  Y'  (t)  -  (y  (“)  •  Y' x'  exists  and  y  (“)  »  0. 

t  *-*■•** 

Proof : 

From  the  arguments  in  the  proof  of  the  last  theorem  we  already  know  the  existence 

and  uniqueness  of  a  solution  on  some  interval  (-“,t^l.  In  order  to  prove  that  the 

solution  exists  globally  in  time,  it  is  more  convenient  to  look  at  (2.2)  rather  than 

at  thi>  equivalent  equation  (2.3).  Solutions  of  (2.2)  continue  to  exist  as  long  as  y 

stay.-  away  from  zero  or  infinity.  From  the  second  and  third  equation  of  (2.2)  one 

obtains  positive  lower  bounds  for  j  q(x)du(x)  and  J  h(X)du(l)  in  every  finite  time 

interval,  provided  that  y  remains  positive,  and  these  bounds  do  not  depend  on  any 

estimate  for  y.  Hence,  if  y  becomes  too  large,  y^-  J  g(X)du(X)  will  dominate 

>ver  fy 1  and  also  over  j  MDdi.(')  (the  latter  beinq  less  than  some  constant  times 

max  y  (  t  ) ) .  Analogously,  if  y  becomes  too  small,  /  h(x)du(X)  will  be  the 
'  <  (-•*>,  tl 
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dominant  term.  It  la  isnadiate  from  this  that  y  cannot  go  to  zero  or  infinity  in 
finite  time,  and  therefore  the  solution  exist  globally. 

For  t  >  tQ,  we  now  have  f  •  0,  and,  putting  a(X)  •  y(X)  -  jj-,  6(M  “  4(1)  - 
we  find  from  (2.1) 


(2.8) 


-  (/  (a  (X)  -  B(X))dy(X)}2 

As  we  know  that  >(1)  and  4(A)  stay  positive,  the  denominators  a  +  and  8  + 
are  always  positive,  and  the  left  side  of  (2.8)  is  therefore  the  derivative  of  a  posi¬ 
tive  function  that  decreases  along  trajectories.  (It  is  easy  to  prove  that  a  and  6 
are  nice  enough  for  all  the  integrals  to  make  sense,  namely,  one  sees  from  (2.2)  that 

»g(M  and  Mi  ( X )  and  hence  Xy(X)  and  X<S(X)  are  bounded).  As  a  consequence,  a 

2 

and  converge  to  zero  exponentially  as  t  -»  •>  in  the  L  -norm  and  a  fortiori  in  the 
L* -norm.  From  the  first  equation  of  (2.3)  one  sees  then  that  y  converges  to  zero 
exponentially,  whence  y  must  converge  to  a  lisdt  exponentially.  Moreover,  one 

easily  concludes  from  the  second  and  third  equations  of  (2.3)  that  a  and  B  in  fact 

«  2 

converge  to  zero  in  the  L  -norm  and  not  only  in  the  L  -norm.  This  concludes  the  proof. 


Remark: 

It  is  almost  trivial  to  prove  ( 7 }  that  y  (•)  > 1  if  f  >  0  and  y  (“)  <1  if 
f  -  0.  Since  the  equation  under  study  describes  the  evolution  of  the  length  of  the 
filament,  if  inertia  are  neglected,  this  Is  a  result  that  one  would  obviously  expect, 
we  have  no  analogue  yet  for  the  full  problem  (0.1),  (0.2). 


-  !  S- 


3.  LOCAL  TIME  EXISTENCE. 


I 


We  now  turn  to  the  study  of  (0.1).  According  to  what  we  have  seen  in  the  last 
chapter,  we  consider  u^ft)  -  b(t)  >  0  as  being  given  at  x  ■  +1,  where  g  is  a 
smooth  function  of  t.  we  want  to  reformulate  (0.1)  in  such  a  way  that  it  fits  into 
the  theory  of  quasilinear  parabolic  equation.  For  this  purpose  we  make  the  following 
substitutions 


p  «  u 


q  »  u 


g.  (X)  -  /  e‘X(t"8)(u  (s)  -  u  (t))ds 

1  1  XX 


,(*>  -  /*  e-x(t-8,pp^ 

lu  <s) 


(s)  u  (t)u  (s) 

XX  X 


U4(t> 

X 


ds 


,(X)  -  / 


-X(t-s)  -2  -2...  .  . 

e  (u  (s!  -  u  (t))ds 

X  X 


g4(X)  »  / 


-X(t-s) 


(u  (s) 

1  XX 


u  (t)u  (t) 

g-g-- — )ds 

u  (s) 

X 


Equation  (0.1)  now  assumes  the  following  form: 


p  =  r  • 
x 


q  *  r 


or  «  3np  2r  -  6gp  3qr  -  2p  /  g  (X)du(X)  — y-  /  g.U)dLi(X) 


g1(X)  =  -Xg^X)  -  ~ 


(3.1) 


g2(X)  -  -Xg2(X)  -  Jf  (g^X)  +E)  ♦  (g^X)  +  E) 

p  p 

2r 

g  (X)  -  -Xg  (X)  +  —5- 
Xp3 

g4(i>  -  -Xg4(X)  -  rxxp2(g3(X)  +  -^y)  -  2rxqp(g3(X)  ♦  -iy) 

Xp  Xp 


with  boundary  condition  p  -  b(t),  rx  «  b(t)  at  x  »  +1. 
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Since  p  »  rx,  the  first  boundary  condition  follows  from  the  second,  once  it  is 

satisfied  initially,  and  we  shall  ignore  it. 

We  will  show  that  (3.1)  can  be  treated  by  the  Sobolevskii  theory.  For  this  we  first 

introduce  some  notations.  H  will  denote  Sobolev  spaces  of  functions  on  (-1,1),  and 
s  Ic  )c 

L  (u,H  )  will  denote  the  space  of  H  -valued  functions  defined  on  C,  which  are 

s-integrable  with  respect  to  the  total  variation  of  p  in  the  Bochner  sense  (for  a 

2  12  s  14 

precise  definition,  see  e.g.  [10]).  We  put  X  »  H  x  (H  )  x  (l  (Uih  ))  .  Moreover, 

s 

in  (3.1)  we  substitute  r  =  r  -  b(t)x  and  introduce  the  abbreviation 
y  -  (p,q,r,gl(g2,g3,g4) .  We  rewrite  (3.1)  in  the  form 


y  =  A(y)y  +  f (y,t)  (3.2) 

where  A(y)  is  defined  as  the  following  linear  operator 


A(y)y '  »  <i\r*  ,  -  (3np_2£’ 

X  XX  p  *  XX 


6np  3grM,  -  Xgpx)  -  ^  F,  -  Xgj  (X)  -  (g1<>) 


r  ' 
£)  xx 
X  4 


P 

+  1|  (g  U,  +E)r;,  -  Xg-  (X)  +J-  i'  -  Xg^  ( X)  -  p2  (,  ( X)  ♦  -±j>  i'  -  2qp  (,  (X) 
p  Xp  Xp  Xp 

with  the  boundary  conditions  r^  *  0  at  x  =  +1. 

We  shall  show  that  (3.2)  satisfies  all  the  requirements  of  the  Sobolevskii  theory 

when  regarded  as  an  evolution  problem  in  ( 1  s  <  °°) .  More  precisely,  we  shall 

prove 

Theorem  3.1: 

Let  1  <_  s  <  <»  be  arbitrary.  Let  =  (p0,q0 , r^, ^)  (  be  given  such  that 

3  -  si 

r.  (  H  ,  r„  =0  at  x  =  +1,  X g .  „  f  L  (u,H  )  and  min  t>  (x)  '•  0.  Then,  for 

°  °'X  -  *'°  xe  (-1,1)  0 

some  T  •>  0,  equation  (3.2)  has  a  unique  solution  y  (  C  ((0,T1,X  )  such  that 


y(0)  -  yQ. 


Proof: 


We  shall  deduce  the  result  from  Theorem  16.2  in  12]  (Theorem  7  in  (BJ  resp.). 
^or  this  we  have  to  verify  the  followinq  conditions  stated  in  (2?  as  (Fl)  and 

(Fl)  The  operator  An  =  A(y  )  is  rtensely  defined,  closed  and  qenorates  an 
analytic  semi q roup. 
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(F3)  For  v,w  in  a  neighborhood  of  yQ  in  %g  there  exists  an  appropriate 
S  t  &  such  that 

||(A(v)  -  A(w))(A(v)  +  6)  1  j|  f.c||v  -  w || 
with  some  constant  c  independent  of  v  and  w. 

(F4)  For  v,w  in  a  neighbourhood  of  y^  and  t,T  e  [0,T]  there  is  some  constant 
c  such  that 


||f  <v,t)  -  f(w,x)  ||  <_  C(|t  -  t|  +  || v  -  w||) 


(FS)  yQ  e  D{Aq) 

(The  conditions  in  [2]  are  more  general,  and  we  have  only  formulated  the  special 
case  applying  to  our  problem) . 

(F4)  and  (F5)  are  trivial  consequences  of  the  smoothness  of  b  and  our  assumptions 

3 

on  the  initial  data.  (F3)  is  clear,  if  it  is  proved  that  the  H  -norm  of  the  r-component 
of  (A(w)  +  6)  can  be  estimated  by  ||y|j.  This  will  be  immediate  from  the  arguments 
leading  to  (FI)  with  yQ  replaced  by  w. 

To  prove  (FI),  consider  the  equation  (AQ  +  B)y  »  y' .  In  the  r-component  this 
leads  to 


3n  -2  -  ,  -3  '  *  * , 

—  p„  r  -  6np_  +  Br  »  r' 

0  0  XX  0  ’0  X 

and  the  equations  for  the  other  components  can  be  trivially  resolved  once  r  is  known. 
It  is  now  a  simple  consequence  of  Theorem  19.2  in  12]  (which  is  due  to  Agmon  and 
Nirenberg  19])  that  if  6  is  in  a  sector  not  containing  the  positive  real  axis,  and 
6  |  is  large  enough,  we  have  an  estimate  of  the  form 

II til  lsl1/2l!t(!  2  *  lei  IIHI  x  <c||r-||  ,  • 

H  H  H  H 


This  concludes  the  proof. 


4.  SOLUTIONS  FOR  SHALL  FORCES. 

The  goal  of  the  present  chapter  is  to  establish  an  analogue  of  Theorem  2.3  for 
the  equation  (0.1),  i.e.  to  prove  existence  of  solutions  globally  in  time  for  small 
forces  f.  f  =  0  now  corresponds  to  the  boundary  condition  ux(t)  ■»  b(t)  =  1  *  b(t)  »  0. 
In  this  case  (3.1)  has  the  trivial  solution  p  =  1,  q  =  0,  r  =  0,  g^  (X)  «  0.  As  a 
first  step  we  shall  study  the  linearization  of  (3.1)  at  this  trivial  solution  with 
homogeneous  boundary  conditions  rx  =  0.  The  linearized  equation  reads  as  follows: 

P  “  rx 
q  -  r 


*  -  T-  -  7  /  g,<*)dp(X)  -  j  /  g4(X)dU(X) 


p  "xx  p  '  *2' 
r 

gx  (X)  =  -Xg1(X)  -  -f- 


g2(X)  -  -Xg2(X)  -  -f- 
2r 

g3  (X)  -  -Xg3(X)  +  ~y 


g4(X)  =  -Xg4(X)  -  -Si 

we  abbreviate  (4.1)  in  the  form  y  =  Ay.  We  shall  study  the  spectral  properties  of 

A  as  an  operator  in  the  space  Xg  of  §3  (1  s  <  “>  is  again  arbitrary) .  Consider 

the  resolvent  equation  (A- a)  y  =  f=  (f  ,f  , f  ,f  (X) , f  (X) , f  (X) , f  (X) ) .  If  -a  is 

.L  2  3  4  j  o  / 

not  in  the  support  of  u,  this  equation  is  immediately  resolved  with  respect  to  p,q 

and  g^(X),  yielding  the  following  equation  for  r 

i  7  f -  (X)  .  f  (X) 

la  r  +  r  .  /  -  1  du(X)  -  ar  =  f  -  -  /  -r— —  dp(X)  -  -  /  r dp(X>  . 

pxxPxx  X(X  +  o)  3  p  ^  X  +  a  pX  +  a 

As  noted  in  §2,  /  -  1 — r-  du(X)  has  a  positive  real  part  for  Re  a  >  0.  Moreover, 

'  X  v  x  +  a)  — 

this  expression  obviously  goes  to  zero  like  yiy  if  a  -*  “  in  any  sector 
(a  e  E  |  -tt  +  ^  +  e  <_  arg  ■f  being  the  angle  of  assumptions  (a)  and  e 

any  positive  number.  From  these  properties  it  can  easily  be  seen  that  the  following 
holds : 
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Proposition  4.1; 


A  is  the  generator  of  an  analytic  semigroup.  Moreover,  the  spectrum  of  A 
consists  of  the  semi-simple  eigenvalue  0  and  a  remainder  lying  strictly  in  the  left 
half  plane. 

Semi-simple  here  means  that  the  resolvent  has  a  simple  pole  at  0,  or  equivalently, 
that  R(A)  #  N(A)  =  Xg,  R(A)  and  N(A)  denoting  the  range  and  nullspace  of  A. 

For  technical  reasons,  the  spaces  X°,  Y°  of  §2  are  not  quite  appropriate  for  the 
study  of  our  present  problem,  and  we  shall  use  the  following  spaces,  which  are  defined 
in  a  very  similar  manner. 

Definition  4.2; 

Let  2  be  a  Banach  space.  Then  Hn(3R,Z)  denotes  the  spaces  of  all  functions 
K  -»  Z  whose  first  n  derivatives  are  square  integrable  in  the  sense  of  Bochner.  Let 
moreover  be 

X°(Z)  =  {v  €  Hn(IR,Z)  |  eatv,e~atv  £  Hn(3R,Z)} 

Y°(Z)  =  {v  :  3R  ->•  Z  I  e"0tv  £  Hn(E,2),3  v  £  Z  such  that  e+0t(v  -  v  >  £  Z)  }  . 

n  1  ®  °° 

Natural  norms  in  X°  and  Y°  are  defined  in  an  analogous  way  as  for  XCT,Y°.  The  use 
n  n  ^  1  n  n 

of  these  definitions  lies  in  the  following  lemma : 

Lemma  4,3: 

Let  the  space  Xj  and  the  operator  A  be  as  above,  and  let  a  >  0  be  small 
enough.  Then  the  operator 

y(t)  -  (a  - 

is  bounded  from  X°(XJ  into  Ya  (N  (A) )  ffi  X°(R(A)  n  D(A) )  ,  where  N  (A) ,  R(A)  and  D  ( A) 
n  2  n  n 

denote  the  nullspace,  range  and  domain  of  A,  resp. 

For  the  proof,  note  that  since  X^  is  a  Hilbert  space,  the  norm  in 
can  easily  be  expressed  in  terms  of  the  Fourier  transform,  thus  reducing  the  statement 

of  the  lemma  to  estimates  on  the  resolvent  of  A.  The  latter  follow  from  Proposition 

O  ■'O 

4.1.  (It  is  this  argument  that  fails,  if  X  is  chosen  rather  than  X  ). 


i 


-20- 


With  these  preliminaries,  it  is  now  easy  to  establish  an  analogue  of  Theorem  2.3 
for  the  nonlinear  problem  (3.1).  Again  we  put  r  «  r  -  b(t)x,  and  we  put 
y  =  (p  -  l,q,r,g1(g2,g3,g4) .  Then  (3.1)  has  the  form 

y  =  Ay  +  f  £y  ,b(t) )  (4.2) 

where  A  is  the  operator  studied  above,  f  is  a  smooth  mapping  from 

ya(N(A))  ©  X°(D(A)  n  R(A))  x  Xn  ( IR)  into  X^tX.)  for  any  n  >  1,  and,  according  to 
n  n  a  n  2  — 

Lerana  4.3,  -  a]  *  is  (for  a  small  enough)  a  bounded  linear  mapping  from  X°(X2) 

into  Y°(N(A))  ®  XJ(D(A)  n  R(A) ) .  The  following  result  is  now  immediate  from  the 
n  n 

implicit  function  theorem: 

Theorem  4.2: 

Let  a  >  0  be  small  enough.  Then,  in  a  neighbourhood  of  y  =  0,  b  =  0  in 
Y°(N(A))  ©  X  (D(A)  n  R(A))  *  Xa(lR),  equation  (4.2)  has  a  unique  resolution  y  =  y(b). 
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ABSTRACT  (continued) 


where  u(x,t)  is  a  real  valued  function  of  x  e  [-1,1]  and  t  e  IR,  with 
the  boundary  condition 


3n 


.  -  iU  (t)  u  (s) 

5  (-i)  *  / 


3t 


u^(s)  u^(t) 


ds  =  f(t) 


at  x  =  +1.  This  equation  is  derived  as  a  model  for  the  elongation  of  thin 
filaments  of  polymeric  liquids,  u  denoting  the  position  of  a  fluid  particle 
in  space,  a  the  memory  kernel,  and  f  the  force  acting  on  the  ends  of  the 
filament.  We  study  the  evolution  of  u,  assuming  the  initial  condition 
u(x,t  =  -°°)  =  x.  It  is  shown  that  under  appropriate  conditions  on  a  and 
f  the  boundary  condition  can  be  uniquely  resolved  with  respect  to  u^.  The 
full  problem  is  transformed  in  such  a  way  that  it  is  approachable  by  the 
Sobolevskii  theory  of  quasilinear  parabolic  equations.  This  yields  the 
existence  of  solutions  to  the  initial  value  problem  on  sufficiently  small 
time  intervals.  Moreover,  we  show  that  if  f (t)  converges  to  zero  exponen¬ 
tially  as  t  -*■  and  is  small  in  an  appropriate  norm,  there  exists  a  solu¬ 
tion  globally  in  time,  which  approaches  a  stationary  limit  as  t  >  +■*. 
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